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i-p^ ' Abstract. - We study the behavior of one dimensional Kac spin glasses as function of the 

f— ! ■ interaction range. We verify by Montecarlo numerical simulations the crossover from local 

Q ' mean field behavior to global paramagnetism. We investigate the behavior of correlations and 

O I find that in the low temperature phase correlations grow at a faster rate then the interaction 

range. We completely characterize the growth of correlations in the vicinity of the mean-field 

^^ ' critical region. 

00 

^^ ' Introduction. - Spin glasses are well understood at the mean field level. After more 

^^ [ then twenty years from the physical clarification of the nature of the spin glass phase of the 
Sherrington-Kirkpatrick model [1], recent progress in mathematical physics [2,3] is rapidly 
leading to a complete mathematical confirmation of the physical implications of replica sym- 
j^ ' metry breaking (RSB). Unfortunately, as soon as one goes beyond mean field, our ability to 
make predictions becomes more limited. Large theoretical efforts devoted to extend the replica 
symmetry breaking theory to finite dimensional systems [4], have not led to an unanimous 

' O I consensus on the nature of the spin glass phase in finite dimension. Approximate renormal- 
ization schemes [5] and phcnomcnological theories commonly known as droplet models [6], 
fj I suggest that glassiness in finite dimension could be very different from the mean field, and 
indeed much simpler. Rigorous attempts to describe low temperature spin glasses are fully 
compatible both with the replica symmetry breaking scenario and with droplet like spin glass 

J^ ' phases [7]. Numerical simulations [8] while giving strong indications that replica symmetry 
H I braking might extend down to dimension three, do not solve the controversy, since it is always 
possible to argue that the systems are not large enough, the samples are not equilibrated etc. 
It has been recently suggested that the nature of the low temperature finite D spin glasses 
can be studied in an asymptotic expansion around mean field using models with long but 
finite interactions of the Kac kind [9]. Some progress has been achieved, proving that for 
large enough interaction ranges, the free-energy is close to the mean field. Moreover, if one 
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admits a mild hypothesis of stochastic stabihty of the Gibbs m.etastate w.r.t. some random 
perturbations, mean-field spin glass order holds at least on a local level: the distribution of 
overlaps on scales of the order of the interaction range is close to its mean field limit. The 
problem about the nature of the spin glass phase in finite dimension can be rephrased as the 
question if the mean-field order can become long-range above a finite lower critical dimension. 
It is clear that in low enough dimension (£' = 1,2) long range order is not possible and local 
mean field order should cross over to paramagnetic behavior on large scales. Even if one feels 
that this case is much simpler that the high dimensional one, a theoretical approach to this 
cross-over is at present still to be developed. In this letter we initiate the study the crossover 
from spin glass to paramagnetic behavior in a ID spin glass model with variable interactions 
via Monte Carlo numerical simulations. 

We first verify in explicit simulations the expected property of genericity of the unperturbed 
model, and show that local overlaps approach mean field behavior. Then we investigate the 
growth of the correlation length for spin-glass order with the interaction range, and we find 
that, analogously to non-disordered models, the correlation length grows more rapidly than 
the interaction range. 

The model. - The model we consider consists in a chain of spins (Ji (i ~ 1, ...,L) with 
periodic boundary conditions, interacting through the Haniiltonian 

M 

H\a\^^^J^ai^aj>. (1) 

where M = {z/2)L, the indexes {i^,j^) are chosen independently from term to term with 
uniform probability among the couples such that |z^ — j^| < -R for some R while the J^ are 
i.i.d.r.v. equal to ±1 with equal probability [10]. If i? = L the model reduces to the Viana-Bray 
diluted spin glass, that for z > 2 admits a low temperature mean field spin glass phase [11]. 
The phase diagram of the model is very simple: for finite R, the model is paramagnetic at 
any positive temperature. On the other hand it has been shown [12] that in the Kac limit, 
R —f oo after the thermodynamic limit L —^ oo, one recovers the mean- field phase diagram, 
with a second order phase transition at a temperature Tc — l/tanh^^{^/z), below which the 
system is in a spin glass phase with full RSB [11]. 

Results. - In order to characterize the behavior of the system we study the local overlap 
between configurations on a scale R. We partition the line {1, ..., L} into disjoint, contiguous 
boxes Bx, {x = \, ...,L/ R) of size R and consider the local overlap between spin configu- 
rations (Ti, and Ti as: qx{<J,T) = iX^ies f^iTi- A simple generalization of the proof given 
in [13] shows that if one couples the original Hamiltonian with suitable finite range random 
perturbations, one generically has that the probability distribution of the local overlap, in- 
duced by the Boltzmann distribution on spin configurations a and r and and the quenched 
random couplings, is close to the overlap distribution function (ODF) of the infinite range 
Viana-Bray model for the same temperature and value of z. This has the characteristic shape 
of mean-field spin glasses with full RSB, with two delta peaks at the extremal values zLqea 
and a smooth part in between [1]. "Generically" refers here to the fact that the property is 
proven almost everywhere in an interval of values of the couplings with the perturbations. 
The stochastic stability property amounts to say that the case of zero couplings is not a sin- 
gular exception. In order to check that indeed this is the case, we simulated the model in ID 
for a value of z = 3, where the mean field critical temperature is T^ — 1.5186. In order to 
equilibrate the system for large samples we used parallel tempering [14]. In this way we could 
reach interaction ranges R = 256 for system sizes of L = 8192 without appreciable finite L 
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Fig. 1 - Distribution of the overlap on scales R for T 
overlap PDF approaches the mean-field distribution. 



: 0.714 and i? = 8, 16, 32, 64. Increasing R the 



effects. All quantities we study are averaged over 100 different samples, and we have checked 
the stability of the average, comparing with the average over 50 samples. 

In figure n we show the behavior of the function of Pniq) for T = 0.714 and various 
values of R. It is apparent that in both cases, increasing the interaction range, the function 
-PrX?) approaches the characteristic mean field shape. This contrast with the behavior of the 
PDF of the global overlap, which for a paramagnet in the thermodynamic limit has a single 
delta peak in zero. We then study the crossover to paramagnetic behavior of the overlap on 
increasing the scale of observation. To this scope, we define overlaps as before, but on boxes 
of size £R and study the distribution as a function of £. In figure |21 this is done for R — 16 
and T ~ 0.714, where we see a clear passage from spin-glass behavior at short lengths to 
paramagnetic behavior at large scale. 

In order to study more quantitatively this crossover, we considered the overlap-overlap 
correlation function 



C{x) 



\<lv<lx+y) 



(2) 



where (•) denotes average over the Boltzmann and quenched coupling distribution. Param- 
agnetic behavior means that this function should tend to zero at large distance, and in one 
dimension one can expect the behavior: 



C{x) = C{0)e^p{-x/O 



(3) 



where C(0) = (g^) is the local average of the overlap square, and ^ is the correlation length 

of the system, which is expected to diverge in the low temperature region T < Tc for R —> oo. 

In this letter we concentrate on the critical region where r = /3 — /3c is small. In that case 

one should cross-over to paramagnetic behavior for finite R to mean field critical behavior for 



i? — > cxD, in which case ^ 



—h== and, for positive r, C(0) 

'\t\ 



T . Analogously to the 



well studied case of non-disordered systems [15], the crossover will be described by scaling 
functions: 



C(0) = (g2)(r,i?)=rVTi?") 



(4) 
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Fig. 2 - Distribution of the overlap on scales IR for T = 0.714 and R = 16, I = 2'' with r = 
1,2,3,4,5,6,7. It is apparent a cross-over from mean-field like two peak behavior to paramagnetic 
Gaussian behavior. 
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where the properties of the functions g and h for large and small argument tR°' should be 
compatible with the expected behavior: namely h should go to positive constants for large 
negative values of x while it should behave as \/x for small argument to cut-off the mean-field 
singularity. Analogously, g should go to a positive constants for large positive values of a;, to 
zero for large negative values and behave as -^ for small x. 

The value of exponent a can be guessed through the observation that RSB effects should 
not affect the critical properties of the system for T > Tc- In that case, one can argue that 
the fluctuations of the order parameter are captured by a cubic field theory [16] which in the 
case of the ID Kac model reads: 



F[q] ^RJdz ({yq{z)f - Tq{zf + \q{zf 



(6) 



where A is a positive constant. Simple scaling analysis predicts then the value a = 2/5. 
This field theory should also in principle suggest the behaviour of the function h{x) for large 
negative arguments, but we did not attempt to follow this route. 

In order to confirm our predictions about the critical exponents and determine the scaling 
functions we simulated the model for z = 3 and various values of R and temperatures. We 
first investigated the behavior of the overlap correlation function at the critical temperature 
Tc, in figure 13 we plot the behavior of C(0) as a function of R showing that the expected 
behavior C(0) ~ l/i?"*/^ is very well respected: we find small corrections to scaling that only 
affect the data points for R = A and 8. In the inset, we plot the scaled function P{q), showing 
that the expected scaling holds for the whole probability distribution. We next investigated 
the behavior of the correlation function. In figure0]we show the collapse of C{x)/C{Q) when 
plotted as a function of x/^ with ^ = (1-2/B.) • A_gain, we verify our scaling form, assuming 
small deviations from scaling. Notice that the assumed exponential form for the correlations 
is very well verified for all a; > 1. The same holds at all the temperatures we looked at. 
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Fig. 3 - C(0) as a function of i? for T = Tc and best fit of the kind v{R) = -^ (l + "I)- In the 
inset, scahng plot of the whole local overlap PDF P{q)v{R) versus ;;j7^- 



We then pass to the task of evaluating the scaling functions. In figureElwe plot, for various 
temperatures, C{0)R^^^ as a function of tFL^/^ which gives the function g{x) = x^g{x). As 
expected, g{x) behaves quadratically for large positive values of x. The behavior at large 
negative values can also be understood, since this is the regime where the effective coupling 



constant in the cubic theory should tend to zero and C(0) 



flVI^I 



leading to g{x) 



Analogously we can understand the behavior of the scaling function h{x) (see fig. EJ: for 
small X, h{x) has the expected square root singularity with a different prefactor above and 
below the critical temperature, as it is usual, while it goes to a constant for large negative 
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Fig. 4 - Collapse of C(a;)/C(0) = /(x/C) for T : 

«!/ = 



Tc, R= 16, 32, 64, 128, 256. The collapse is obtained 



for C^ 
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Fig. 5 - Scaling of the C(0, T, R) close to Tc. In the high negative argument region we show the fit 
g(x) = 0.67/\/— a; in the region x < —0.5, which is indistinguishable from the data. 



values. The behavior of the function h for large positive values should describe the cross-over 
from critical to low temperature behavior, in particular, it should determine the behavior of 
the correlation length in the low temperature phase. Unfortunately, the data we have, though 
indicate that a power law behavior ^ ~ R^ may persist at low temperature with an exponent 
CO larger then 1/5, do not allow a precise determination of the exponent lo which would require 
larger interactions ranges R. 

Summary. - Summarizing we have studied the crossover from paramagnetic to mean field 
behavior in a ID Kac spin-glass model. Our work has a qualitative aspect, from which we get 
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Fig. 6 - Critical scaling of the correlation length. The curves are fitted by h{x) — 0.3^, ^ 277-2^ 
negative x and by h{x) — 0.42-^3; — 0.35a:: for positive x. 
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evidence that stochastic stability holds at low temperature, and a quantitative aspect where 
we study the crossover from paramagnetic to mean- field behavior close to Tc- We characterize 
this crossover through scaling functions describing the behavior of the local Edwards- Anderson 
parameter C(0) and the correlation length ^. We find that while in the high temperature phase 
the correlation length is, in units of the interaction range R, independent of R, in the low 
temperature phase a dependence on R sets in. The correlation length grows as R^'^ in units 
of R at the critical point, while it grows faster at lower temperatures. This means that in 
units of lattice constant the correlation length grows as i?^+^/''' or faster. This result shows 
that the rigorous analysis of [13] just provides a lower bound to the size of the regions where 
local mean- field order holds. Further studies will be necessary to make a quantitative analysis 
deep in the low temperature region and to go to higher dimension. 
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